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ABSTRACT 



Context. We re-analyze the precision radial velocity (RV) data of HD 208487, HD 190360, HD 188015 and HD 114729. These stars are 
supposed to host Jovian companions in long-period orbits. 

Aims. We test a hypothesis that the residuals of the 1 -planet model of the RV or an irregular scatter of the measurements about the synthetic 
RV curve may be explained by the existence of additional planets in short-period orbits. 

Methods. We perform a global search for the best fits in the orbital parameters space with genetic algorithms and simplex method. This makes 
it possible to verify and extend the results obtained with an application of commonly used FFT-based periodogram analysis for identifying the 
leading periods. 

Results. Our analysis confirms the presence of a periodic component in the RV of HD 190360 which may correspond to a hot-Neptune planet. 
We found four new cases when the 2-planet model yields significantly better fits to the RV data than the best 1-planet solutions. If the periodic 
variability of the residuals of single-planet fits has indeed a planetary origin then hot-Neptune planets may exist in these extrasolar systems. 
We estimate their orbital periods in the range of 7-20 d and minimal masses about of 20 masses od the Earth. 

Key words, extrasolar planets— Doppler technique— HD 208487— HD 190360— HD 188015— HD 114729— HD 147513 



1. Introduction 

The precision radial velocity (RV) surveys for extrasolar plan- 
ets approach the limit of detection of Neptune-like compan- 
ions. In the past year, 14-Earth mass planet has been an- 
■ nounced bv ISantos et alJ (|2p04). A lo w-mass companion ha s 
been found in the Cancri system ( Mc Arthur et al. "2OO4I. 
. Almost at the same time, a Neptu ne-like planet h as been also 
' discov ered bv lButler et alJ (l2004h . Very recentlv. [Rivera et al] 
(|20P5) announced a discovery of ~ 7 Earth-mass planet in 
~ 2d orbit about t he famous Gliese 876 planetary system. Also 
IVogt et alJ ( l2005t) announced a new multi-planet system about 
HD 190360 involving a Neptune-mass planet in ^ 17. Id or- 
bit with a Jupiter-like companion with ~ 3000d orbital period. 
These discoveries change the view of multi-planet systems. 
The found low-mass planets likely are rocky bodies rather than 
massive gaseous hot-Jupiters. The architecture of their environ- 
ment, involving also giant planets in relatively distant orbits, 
begin to resemble our own Solar system. 

There are already speculations as to which of the extrasolar 
systems involving Jupiter-like planets can contain smaller com- 
panions. At present, such objects cannot be clearly detected be- 
cause of small amplitude of the contributed RV signal and not 
frequent enough observations. Nevertheless, the existence of 
such planets seems to be natural. The dynamical relaxation of 



planetary systems jAdams & Laughliij|2003l) can lead to final 
configurations with orbits filling up the dynamically available 
space of stable orbits (Laskar 2000). In particular, many extra- 
solar systems containing distant Jupiter-like planets may con- 
tain smaller, still undetectable companions, which can e asily 
survi ve in a close neighborhood of the parent star. Baraffe et alJ 
( I2OOS) show that the three recently detected hot-Neptune plan- 
ets (about GJ 436, p' Cnc, yU Arae) may originate from more 
massive gas giants which have undergone significant evapo- 
ration. This work suggests that hot-Neptunes and hot-Jupiters 
may share the same origin and evolution history and they can 
be as frequent as are hot-Jupiters. The numerical simulations 
by iBrunini & Cioncoi ( l2005l) show that the formation of com- 
pact systems of Neptune-like planets close to solar type stars 
could be a common by-product of planetary accretion. In turn, 
these theoretical works are supported by several n ew detections 
of hot-Neptune and hot-Saturn candidates (e.g.. 'Bonfils e t all 
^005; Udrv et al. 2005; Curto et al. 2005; Fischer et al. 2005^^ 
The detection of Neptune-mass exoplanets is a break- 
through in the field but the inevitable problem still lies in the 
limitations of the Doppler technique. The radial velocity sig- 
nal of small, Neptune-like planets revolving close to the host 
star has the semi-amplitude of a few meters per second, i.e., at 
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a similar range as the measurement errors and chromospheric 
variability of the RV. This can be a source of erroneous inter- 
pretation of the observations. For instance, in a recent work, 
[wisdom (2005) clair ns that the orbital period of 2.38d of com- 
panion d detected bv lMcArthur et all ( l2004l) in Cancri sys- 
tem is an alias of two periods originating from the orbital mo- 
tion of planets b and c and, instead, he found an evidence of 
an another planet in ~ 261 d orbit. Actually, a proper inter- 
pretation of the observations is still very difficult. Often, they 
can be equally well modeled by qualitatively different config- 
urations and even the number of planets cannot be resolved 
without doubt. An example will be given in this paper too. 

In this work, for a reference and a test of our approach, 
we re-analyze the p ubl ished RV data of two stars: HD 208487 
llTinnev et alJl2005l) — Ig regor vl ll2005a') claims that a second 
planet acco mpanies a confir med Jovian companion b, and 
HD 190360 IVogt et al.ll2005h — the discovery team claims the 
existence of a hot-Neptune and a Jovian planet in long-period 
orbit. Next, using the same metho d, we study oth er three extra- 
solar systems: HP 188015 jMarcv et alJi2005bl). HP 1 14 729 
llButler etalJl2003l) and HP 147513 ("Mavor et alJ I2004 '). in 
which Jupiter-like companions have been discovered in long- 
period orbits. However, the residuals of the 1-Keplerian fits 
are either relatively large or we see "suspicious" scatter of the 
residuals. We use the available measurements to test a hypoth- 
esis that additional short-period planets exist in these systems 
and their signals may be already hidden in the existing data, but 
likely at the detection limit. Our reasoning follows the theoret- 
ical predictions that Neptune-like planets can easily form and 
survive close to the parent star. The main difficulty is a small 
number of measurements, typically about of 30^0 data points 
spread over several years. It makes it very difficult to determine 
orbits of such planetary candidates without a doubt. In partic- 
ular, this concerns the identification of the orbital periods with 
the application of somehow traditional Fourier-based analysis. 
Instead, we prefer a global search for the best-fit 2-planet con- 
figurations using the genetic algorithms and scanning (^y)'^^ in 
the multidimensional space of orbital parameters. This could be 
an alternative to the already classic Lomb-Scargle periodogram 
(e.g., Butler et al. 2004) or sophistic ated Bayesian-b ased anal- 
ysis, which is recently proposed bv taregorvl l l2005 a b). That 
author claims an evidence of additional planets using the same 
RV measurements which are modeled by 1 -planet configura- 
tions by the discovery teams, about HP 73526 and HP 208487. 

2. The fitting algorithm, numerical setup and tests 

In all cases studied in this paper, we assume that the second, 
putative planet is closer to the star than the already known com- 
panion b in relatively long-period orbit. To model the RV sig- 
nal, we use the standard formulae (e.g., Sm arLl949.) . For every 
planet, the contribution to the reflex motion of the star at time t 
is the following: 

y^{t) = /:[cos(w + v(0) + e COS + Vo, (1) 

where K is the semi-amplitude, to is the argument of pericenter, 
v(f) is the true anomaly involving implicit dependence on the 
orbital period P and the time of periastron passage Tp, e is the 



eccentricity, Vq is the velocity offset. There are some arguments 
that is best to interpret the derived fit parameters {K, P, e, a>, Tp) 
in terms of Keplerian elements and minimal masses related 
to Jac obi coordinates (Lee & Peak 2003; Gozdziewski et al., 
l2003l) . 

To search for the best-fit solutions, we applied a kind of 
hybrid optimization. A single program run starts the genetic al- 
gorithm s (OAs) optimization . In particular, we use the PIKAIA 
code by 'CharbonneauJ (Il995l) . The GAs has many advantages 
over more popular gradient-ty pe algorithm s, for instance, the 
Levenberg-Marquardt scheme dPress et alJ[l992l) . The power 
of GAs lies in their global nature, the requirement of knowing 
only the (xl)^^^ function and an ease of constrained optimiza- 
tion. In particular, GAs per mit to define parameter bounds or 
add a penalty term to O^v)'^^ JGozdziewski et al.l2003[l2005h . 

In our code, the initial population of randomly chosen 
2048 members represents potential solutions to the RV model. 
The code produce 128 new off-spring generations. The best 
fit found by PIKA IA is refine d Iry the simplex algorithm of 
Melder and Nead dPress et alJll992h . In that method, we set 
the fractional tolerance for (yJ)'^^ variability to 10"'^. Such 
procedure runs thousands of times for every analyzed data set. 
Thanks to non-gradient character of both algorithms we may 
bound the space of examined parameters, according to specific 
requirements. Usually, due to a small number of measurements, 
the best-fit is barely constrained and tends to be large, in 
particular in the range of very short periods. Nevertheless, we 
assumed that the tidal circularization limits the eccentricity of 
the short-period orbit and we arbitrarily set its upper bound to 
0.3. The orbital periods of putative planets are searched in the 
range of [2,136] d. 

The internal errors of the RV data are rescaled according 
to cr^ - o"^ + crj, where cr is the joint uncertainty, cr,^ and 
cTj is the internal error and adopted variance of stellar jitter, 
respectively. Typ ically, \ye cho ose crj following estimates for 
Sun-like dwarfs ilWrighlll2005 ) or we set its value as quoted 
by the discovery teams. The particular estimates are given in 
Table 1. 

The best fits found in the entire hybrid search are illustrated 
by projections onto planes of particular parameters of the RV 
model, Eq. 1 . Usually, we choose the (P, K)- and {P, e)-planes. 
Marking the elements within the formal Icr, 2cr, 3cr confidence 
intervals of the best-fit solution (which is always marked in the 
maps by two crossing lines) we have a convenient way of vi- 
sualizing the shape of the local minima of (y?)'^^ and to obtain 
reaso nable estimates of the fit errors ( Bevington & Robin son! 
l2003l) . 

Still, even we find very good best-fit configurations, a prob- 
lem persists: is the obtained fit not an artifact caused by specific 
time-distribution of measurements or the stellar jitter? In order 
to answer for this question we apply the test of scrambled ve- 
locities (Butler et al. 2004). We remove the synthetic RV con- 
tribution of the firstly confirmed planet from the RV data, ac- 
cording to its best-fit parameters. Next, the residuals are ran- 
domly scrambled, keeping the exact moments of observations, 
and we search for the best-fit elements of the second (usually, 
inner) planet. For uncorrected residual signal we should get 
close to the Gaussian distribution of (xlY^^- Forevery synthetic 
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data set, the minimization is performed with the hybrid method. 
Let us note, that in this case we deal with 1 -planet model (and 
6 parameters only) thus the large population in the GA step and 
precise simplex code assure that the best-fit solution can be re- 
liably found in one single run of the hybrid code. The runs of 
the code are repeated at least 20,000 times for every analyzed 
data set. If the best fit corresponding to the real (not scrambled) 
data lies on the edge of this distribution, such experiment assure 
us that the real data produce much smaller (xlY^^ (and an rms) 
than the most frequent outcomes obtained for the scrambled 
data. Such a test provides a reliable indication that the residuals 
are not randomly distributed. A likelihood pu that the residuals 
are only a white noise can be measured by the ratio of the best 
fits with (xiy^^ less or equal to that one corresponding to the 
real data against all possible solutions. 

Unfortunately, even with a positive outcome of the test 
of scrambled residuals we still are not done. Even if we find 
a strong indication that the residuals are not randomly dis- 
tributed, and we find a periodic component in the residual sig- 
nal, then we still cannot be sure that it has a planetary origin 
and it does not eliminate jitter as a source of the RV signal. 
We are warned by the referee (and here we literally follow his 
arguments) that jitter can have some spectral power usually as- 
sociated with the rotation period of the star, but other peaks can 
be present. The spectral structure of very low-amplitude jitter is 
not well-known yet, but for large amplitude jitter in active star, 
the power spectrum can s how several peaks and aliases. One 
example is HD 128311 in lButler et alJ ll2003h : an another one 
is, possibly, HD 147513 jMavor et alJl2004 analvzed in this 
work. Therefore, the noisy periodic signal may be due to jitter 
and it may be possible that some secondary periods, which are 
present in the residual signals, are corresponding to aliases of 
the stellar rotation period reinforced by the noise. This is the 
main reason why our work cannot be considered in terms of a 
discovery paper; additional measurements and analysis are re- 
quired to confirm (or withdraw) the hypothesis about the plan- 
etary nature of the RV variability. Our basic goal is to test the 
approach on a few non-trivial and, possibly, representative ex- 
amples. We try to find characteristics of the RV da ta sets which 
can co nfirm the periodogram-based approach (e.g. JButler et alJ 
l2004l) and, possibly, may be helpful to extend its results. Our 
analysis may be useful to speed-up the detections of new plan- 
ets and planning the strategy of future observations. 

2.1. Test case I: HD 208487 

The first test case of HD 208487 jTinnev etalJ l2005h illus- 
trates nuances of modeling the RV data. The observational win- 
dow of this star spans about of 2500 d but the number of data 
points is only 3 1 . The discovery team found a planetary com- 
panion with minimal mass of 0.7 nij and orbital period about 
of 130 d. An rms of this 1 -planet solution is ~ 7 m/s, while 
the joint uncertainty is about of 6 m/s, assuming jitter estimate 
of 3 m/s, a rea sonable value for the quiet, 6 Gyr old star with 
R'fjj^ ^ -4.90 ("Tinnev et al." 2005). The data have been reana- 
lyzed bv iGregorv. C2005 a). He applied his Bayes-periodogram 
method and found much better 2-planet solution, explaining the 
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Fig. 3. Stability analysis of the orbital fit with putative 
Neptune-like body in the HD 208487 system. The upper panel 
is for the eccentricities, the lower panel is for MEGNO com- 
puted for 3-10^ orbital periods of the outer companion. 

residual signal of the 1 -planet solution by the presence of a dis- 
tant, second companion in 998 d orbit. This solution has signif- 
icantly better rms of 4.2 m/s and (xlY^^ - 0.83 than the best 
single planet solution, (xlY^^ - 1.21. 

The analysis of the HD 208487 data is a good refer- 
ence for our simpler approach. First we did the hybrid search 
for the best 2-planet configurations, assuming that 130 d < 
Pc < 1200 d. The results of this search are illustrated in 
Fig. [0 Curiously, the best-fit has (xl)^^^ - 0.62 and an rms 

3.1 m/s, so i t is apparently better than the fit quoted by 
iGregorvl (l2005al) . The best-fit parameters {K,P,e,(jj,T^ - Tq) 
are: (17.532 m/s, 129.224 d, 0.195, 188°.232, 1421.112 d) and 
(13.117 m/s, 997.953 d, 0.433, 190°.470, 2358.258 d) for the 
inner and outer planet, respectively (where ro=JD 2,450,000) 
and Vo = 9.64 m/s. Note much larger eccentricity of the outer 
planet than in the fit bv lGreg orv (2005a), - 0.2. The Icr scat- 
ter of the derived solutions around the best-fit (Fig.Q provides 
the error estimate of the orbital elements. In general, the pa- 
rameters of the inner planet are well fixed, while the elements 
of the outer body have much larger formal uncertainties. The 
errors of K^-, P^, are about of 100 d, 0.3, and 10 m/s, respec- 
tively. 

Next, we explored the inner space to the ~ 130 d orbit, 
assuming 2 d < < 130 d. Surprisingly, we found ever 
better fit, given in Table 1, with (xl)^^^ - 0.581 and an rms 
^ 2.82 m/s. The orbital elements of the best-fit solutions found 
in the entire search are projected onto the parameter planes and 
they are illustrated in Fig. |2] Note a small scatter of the el- 
ements of the outer planet and two sharp minima in P^. The 
best-fit has f c - 14.5 d and Sc - 0.42. Such large eccentric- 
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Fig. 1. The parameters of the best-fit solutions to the 2-planet model of the HD 208487 RV data projected onto the (P, e)- and 
(f , /T) -plane. In this search, a more distant companion to the formerly known Jupiter-like body in ^ 130 d orbit is assumed. 
Upper plots are for the inner planet, bottom plots are for the outer companion. The values of (xl)^^^ of the best-fit solutions 
are marked by the size of symbols. Largest circle is for (xl)^^^ equal to 0.619 ; smaller symbols are for Icr solutions with 
Cr2)i/2 g [0.619, 0.710), 20- solutions with (xl^'^ e [0.710, 0.841) and 3cr solutions with (xl)^'^ e [0.841, 1.0) (smallest, filled 
circles), respectively. 
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Fig. 2. The parameters of the best-fit solutions to the 2-planet model of the HD 208487 RV data projected onto (P, e)- and 
(P, A')-plane. In this search, an inner companion to the formerly known Jovian planet in ^ 130 d orbit is assumed. Upper plots 
are for the outer planet, bottom plots are for the inner companion. The values of (xlY^^ of the best-fit solutions are marked 
by the size of symbols. Largest circle is for Ctf^)'^^ = 0.581 (the best-fit solution. Table 1); smaller symbols are for solutions 
within Icr interval of the best-fit with Od)^'^ e [0.581, 0.677), 2a- solutions with Od^'^ e [0.677, 0.813) and Sa solutions with 
(xiy^ e [0.813, 0.98) (smallest, filled circles), respectively. 
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Fig. 4. Lomb-Scargle periodogram of the residuals of the 
Jupiter-like companion's with Pb ~ 130 d about HD 208487. 



ity may contradict the requirement of a tidally circularized or- 
bit. However, the orbital evolution of the whole system reveals 
large-amplitude variations of (Fig.|3} which are forced by the 
mutual interactions between companions. We also face an an- 
other trouble, a small periastron-apoastron distance of the two 
putative planets ~ 0.2 AU. Nevertheless, the system appears to 
be dynamically stable. Th e evolution of the MEGNO indicator 
( iGozdziewski et alJ200li) computed over 3 ■ 10^ Pb is illustrated 
in the lower panel of Fig. |3l Its oscillations about of 2 mean 
a quasi-periodic configuration and strictly confirm apparently 
regular behavior of e^.c- 

For a reference, we computed the Lomb-Scargle (LS) peri- 
odogram (Fig. 13 of the residuals of the dominant 1-planet sig- 
nal having the period of 130 d. The largest peaks about of 28 d 
and 998 d (the period of a hypothetic long-period planet) most 
possibly are aliasing with the best-fit period ^ 14.5 d. However, 
this period is favored by the smallest (Yv)'^^- The aliasing of the 
two dominant periods is illustrated in Fig.|5] This figure shows 
the synthetic RV signals of the putative configurations with the 
short- and long-period companions of planet b. Curiously, the 
orbital periods of the two planets in the model with the short- 
period orbit are very close to the 9: 1 ratio. This again could be 
an effect of aliasing, nevertheless the best-fit solution with sin- 
gle planet has P\, about of 14.5 d and much larger rms ~ 7 m/s 
than in the best 2-planet solution. 

Some more interesting results brings the test of scrambled 
residuals, Fig.|6l It appears that the real residual signal lies far 
from the maximum of the Gaussian-like distribution and it is 
unlikely that the residuals are only a white noise (pu ^ 3 ■ 10"^). 
Also Kb ^ 10.5 m/s is much larger than the joint uncertainty of 
the measurements. Because HD 208487 is a very old and quiet 
G-dwarf with low chromospheric activity jTinnev et al.k2005.) . 
the hypothesis of a hot sub-Neptune planet in this extrasolar 
system seems to be very appealing. A peculiarity of this system 
could be the proximity to the 9: 1 mean motion resonance of its 
companions. 
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Fig. 6. A histogram of Otv)'^^ for Keplerian fits to 33,000 sets 
of scrambled RV residuals of the synthetic signal of the Jupiter- 
like planet in ~ 130 d orbit (HD 208487, see Table 1). The 
position of the best-fit to the real data is marked with an arrow. 
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Fig. 8. Lomb-Scargle periodogram of the residuals of the long- 
period signal which is seen in the HD 190360 data. 

2.2. Test case II: HD 190360 

In a recent paper, IVoet et all l Eool published the precision RV 
data for HD 190360. The variability of the RV is explained by 
the existence of two planetary companions: a Jupiter-like planet 
in long-period orbit, Pb ~ 3000 d, and a Neptune-like body 
in ^ 17.1 d orbit. The authors conclude that the existence of 
the smaller planet is the best explanation of the residual signal 
of the single planet solution. Its contribution has small semi- 
amplitude K ^ 4.6 m/s. Yet the orbital period is about of a 
half of the rotational period of the star, 36-44 d, and the RV 
variability may be an effect of spot complexes on the stellar 
surface ( Vogt et al. 2005). 

The elements of a few hundred fits gathered by the hybrid 
search are illustrated in Fig. (the elements of the best-fit are 
given in Table 1). The period of the inner planet corresponds to 
a single, very deep minimum of {xiy^^- It is clear, that has a 
lar ge error ~ -15. Our fit is slightly better than that one quoted 
bylyogt et al,i ( 2005) but - 0.2 is also different as well as 
we have got slightly larger ^5.2 m/s. Note that we obtain a 
similar uncertainties of 100 d for Kb and of 0.05 for e^- 

The Lomb-Scargle periodogram (Fig. |8} of the residuals 
of the long-period orbit reveals a strong, dominant peak at 
17.1 d. This an excellent confirmation of the results of the hy- 
brid search. The test of scrambled residuals reveals that the 
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Fig. 5. Synthetic RV signals of the two solutions to the RV data of HD 208487. The thick line is for the outer companion to the 
known 130 d planet. The thin line is for the putative inner companion. The original RV data are marked with error-bars. For a 
clear illustration, only a part of the second curve is plotted. 
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Fig. 7. The parameters of the best-fit solutions to the 2-planet model of the HD 190360 RV data projected onto the (P,e)- and 
(f , ^r)-plane. Upper plots are for the outer planet, bottom plots are for the inner companion. The values of (xl)^^^ of the best- 
fit solutions are marked by the size of symbols. Largest circle is for (xl)^^^ - 0.679 (the best-fit solution. Table 1); smaller 
symbols are for Icr solutions with C^v)'^^ e [0.679, 0.70), 2cr solutions with (xlY'^ e [0.70, 0.74) and 3cr solutions with (xl)^'^ e 
[0.74, 0.79) (smallest, filled circles), respectively. 



residual signal is hardly random. The maximum of {xiy^^ for 
scrambled data is shifted by ^ 0.4 from the value of (xlY^^ for 
the real data. 

Now, having so clear picture, we analyze a smaller sample 
of measurements. To perform this test, we selected 38 first mea- 
surements spanning ~ 2200 d. The number of data points and 
their time-coverage are similar to other cases which we study 
in this work. The described above procedure was repeated on 
such synthetic data set. 



Still, the best-fit period of the inner planet can be clearly 
identified (Fig.llO> but the semi-amplitude Ac has been slightly 
shifted. The elements of the outer planet are determined with 
much larger errors than the uncertainties obtained for the full 
data. In particular, the error of is about of 350 d. The domi- 
nant period is still well recognizable in the Lomb-Scargle peri- 
odogram, Fig.^J Finally, the test of scrambled residuals also 
give a clear picture, albeit the histogram (Fig. I12> has much 
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Fig. 10. The parameters of the best-fit solutions to the 2-planet model of the RV data of HD 190360 limited to the first 38 
measurements, projected onto (P, e)- and {P, /r)-plane. Upper plots are for the outer planet, bottom plots are for the inner com- 
panion. The values of (xl)^^^ of the best-fit solutions are marked by the size of symbols. Largest circle is for (xl)^^^ - 0.733; 
smaller symbols are for Icr solutions with (^2)i/2 g [0.733,0.793), 2cr solutions with (xlY'^ e [0.793,0.89 and 3o- solutions 
with (xiy^^ e [0.89, 1.00) (smallest, filled circles), respectively. Vertical and horizontal lines mark the best-fit obtained for the 
full data (see Fig.0and Table 1). 
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Fig. 9. A histogram of ixlY^^ for Keplerian fits to 20,000 sets 
of scrambled RV residuals of the synthetic signal of the outer 
planet in HD 190360 (see Table 1). The position of the best-fit 
to the real data is marked with an arrow. 
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Fig. 11. Lomb-Scargle periodogram of the residuals of the 
long-period signal in HD 190360 when only 38 first measure- 
ments are analyzed. Compare with Fig.|8lfor the full data set. 



larger width and (xl)^^^ of the real data is closer to the maxi- 
mum of the distribution. 

Even for a half of the available measurements, the peri- 
odic nature of the residual signal of the single-planet fit can be 
clearly recognized. One might conclude that the inner planet 
around HD 190360 might be detected about 2-3 years ago. 

The results obtained for the HD 190360 data provide a valu- 
able test of our method and a good reference point to further 
analysis. Let us recall, that the model with two planets yields 
low residuals with an rms of 3.5 m/s and ixlY^^ much smaller 



than 1, similarly to the three systems which are described be- 
low. 

3. HD 188015 

Recently, iMarcv et al.l (l2005bl) published the precision RV of 
HD 188015. The variability of the RV is explained by the pres- 
ence of a planetary companion with minimal mass 1.26 mj 
and orbital period ~ 1 .25 yr Our attention has been drawn on 
large rms of this 1-planet solution, ^ 8.6 m/s. The mean error 
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Fig. 12. A histogram of (xlY'^ for Keplerian fits to ~ 27, 000 
sets of scrambled RV residuals of the synthetic signal of the 
outer planet in HD 190360 for the first 38 measurements. The 
position of the best-fit for the real data is marked with an arrow. 
Compare with Fig.|9l 
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Fig. 15. Lomb-Scargle periodogram of the residuals of the 
long-period signal visible in the RV data of HD 188015. 



of the RV data is about of 4.1 m/s and expected stellar jitter 
^ 4 m/s (Marcv et al. 2005b). The mean measurement error 
added in quadrature to the jitter is about of 5 m/s, still much 
smaller than the rms of 1 -planet solution. However, this can be 
partially explained by a linear trend apparently present in the 
RV signal. If one adds such a trend to th e fit model, the rms 
goes down by 1 m/s ( Marcv et a l. 2005b, also our Fig. I13h). 
The trends in the RV data frequently indicate a distant body 
(e.g..lMarcv et alJl2005atlJones et alJl2002HGozdziewski et all 
l2003 >). Nevertheless, we assumed that the second putative body 
is not a distant companion, but rather a planet closer to the star 
than already detected companion b. 

The elements of the best-fit solution found in the entire hy- 
brid search are illustrated in Fig. (panels b,c,d) and its el- 
ements are given in Table 1. The synthetic Doppler signal is 
plotted over the data points in Fig. 113b . Surprisingly, the ap- 
parent slope of the RV data can be only an artifact related 
to a specific time-distribution of measurements. The value of 
(xiy^^ - 1-021 indicates a statistically good solution. Its rms 
is about of 5 m/s, consistent with the adopted estimate of the 
mean uncertainty. The semi-amplitude of the inner planets' sig- 
nal is ~ 7 m/s. Assuming the mass of the host star ^ 1.08 Mq 
(iMarcv et al. 2005b), we derived its minimal mass about of 
0.07 OTj and semi-major axis about of 0.08 AU. Thus, the new 
putative object can be classified as a hot-Neptune planet. 

According to the Icr-level of (^^)'^^, marked in Fig.^^ the 
uncertainty of e\, is ~ 0.1, and of A'b is ~ 5 m/s. Similarly, 
the Icr uncertainties of the elements of the inner planet can be 
estimated by ~ 0.3 and 1 m/s, respectively. There are visible 
several local minima in the (P\,, eb)-plane, which are very nar- 
row with respect to P^, about of 0.05 d at the Icr level. It is clear 
that the eccentricity of the inner companion (also its argument 
of periastron — not shown here) cannot be well fixed. 

A local minimum at P^. ^ 15.37 d which also well mod- 
els the observations (at the Icr confidence interval of the best- 
fit solution) may correspond to an alias of the best-fit period 
8.386 d. It could be also an alias of the rotational period ~ 30 d 
of the star, assuming that it is smaller than the 36 d estimate 
by (Marcy et al. 2005b) . Indeed, there is a minimum of (xl)^^^ 



at the Icr-level of the best-fit about of 30 d and ~ 0.3, see 
Fig. HI 

The results of the hybrid search perfectly agree with the 
LS periodogram (Fig. I15> of the signal after removal of the 
contribution from the outer planet. All the relevant periods 
can be found in the periodogram too. This assures us that the 
GA search is robust in locating the best-fit solution even in 
the case when O^v)'^^ has several equally deep local minima. 
Additionally, we did an independent gradient-based scan over 
f c and that test also confirmed the excellent performance of the 
GAs. 

The test of scrambled velocities which is illustrated in 
Fig.Elhelps us to explain the trouble of finding a clear, sec- 
ondly dominant short-period in the data. The histogram of 
(xiy^^ has the maximum about of 1.1 and (xl)^^^ - 0.93 of the 
real residuals is not well separated from this maximum. The es- 
timated Pa - 0.03 is significant and indicates that the residual 
signal is rather noisy. Yet the best 2-planet fit improves signifi- 
cantly the 1 -planet solution. 

We are warned by the referee (and we quote him again here) 
that the jitter estimate by Marcv et al. (2005b) comes from the 
study of jitter levels in the Carnegie planet search sample by 
IWrightl il2005h . In that article, it turns out that at the activity 
level of HD 188015 corresponds to an average jitter of the 
quoted 4 m/s, but with a wide dispersion from case to case, 
and that levels of 6-7 m/s are not at all uncommon. In fact, 
for low values o f the C alcium index there is little activity-jitter 
correlation in Wright (2005), and any values between 1 and 
10 m/s are almost equally likely. This may invalidate our ar- 
gument about the excess scatter and the value of (xlY^^ in the 
2-planet model. Nevertheless, what we can do is to adopt an 
uniform jitter variance for all measurements. Then a larger jit- 
ter will only "flatten" (xlY^^, keeping its minimum at similar 
position in the parameters space. The results for HD 190360 
may be used as an argument supporting the planetary hypothe- 
sis — in that case the 2-planet fit also has an rms about of 1 m/s 
smaller than the joint uncertainty. Yet we encountered a similar 
problem when analyzing the next system, around HD 1 14729, 
described below. 
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Fig. 13. The best-fit solutions to the RV data published by iMarcv et"al] ( l2005bl) for HD 188015. a) 1 -Kepler model with a 
linear trend which is apparently present in the data, b) 2-Kepler model is consistent with the measurement error, ixlY^'^ - 1-02, 
rms^ 4.9 m/s. c) Period-phased RV signal of the inner companion, d) A close-up of the synthetic signal of the best-fit 2-planet 
configuration (see Table 1). 
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Fig. 16. A histogram of (xl)^l^ for Keplerian fits to 20,000 sets 
of scrambled RV residuals of the synthetic signal of the outer 
planet in HD 188015 (see Table 1). The position of the best-fit 
to the real data is marked with an arrow. 



4. HD 114729 

Similarly to most of the stars monitored in precision RV sur- 
veys, the HD 1 14729 is a Sun-like star of G3 V spectral type 
(GOV by Hipparcos). Its RV data indicate a Jupiter-like com- 
panion with minimal mass 0.84 otj. orbi tal period ^ 3.1 yr 
and eccentricity ^ 0.32 dButler et alJ2003l) . An rms of the sin- 
gle Keplerian fit to the data is 5.34 m/s. It is larger than the 
mean of internal e rrors, ^ 3.6 m/s. The expected jitter is 4 m/s 
dButler et al.1 12003). The mean measurement error scaled in 
quadrature with the jitter variance is comparable with the rms. 
Nevertheless, having in mind a significant scatter of the mea- 
surements about some parts of the synthetic 1 -planet RV curve 
(see Fig.ll7h). and the presence of the distant giant planet, we 
assumed (as in the previous case) that there exist an another 



one, inner planet in the system in short-period orbit, [2,136] d. 
To test this hypothesis, we did the same analysis of the RV data 
as described in previous sections. The results are illustrated in 
Fig.llTh.b.c.d which is for the graphical illustration of the best- 
fit parameters. Figure[^is for the Lomb-Scargle periodogram 
of the residual signal of the long -period companion. Figure [T9l 
illustrates the shape of (xl)^^^ about the best-fit and the distri- 
bution of the best fits obtained in the hybrid search. The or- 
bital parameters of the best fit are given in Table 1. Clearly, 
the results of the hybrid search agree perfectly with the LS 
periodogram — the dominant peak about of 7 d is clearly vis- 
ible in the (f c, ec)-plane as the deepest minimum of (xl)^^^ in 
the region of short periods. The best-fit eccentricity of the in- 
ner planet is found close to the assumed limit of 0.3. However, 
in a few tests with smaller jitter (~ 3 m/s) we obtained better 
bounded minimum of (xl)^^^ and smaller best-fit - 0.2, so 
this element is barely constrained by the RV data. 

An inspection of the LS periodogram shows some other rel- 
atively strong peaks, for instance about of 180 d, which most 
likely are aliasing with the period of 13.844 d. 

FigureEUlis for the histogram of (xl)^^^ obtained for scram- 
bled residuals of the outer planet's signal. The probability that 
the residual signal of the outer companion has a random origin 
is very small, Ph - 0.003. The 0^^2)1/2 foj. t^e best-fit solution 
has well defined minimum about of 0.77. Its rms of ~ 3.6 m/s is 
very close to the mean of the internal errors. Assuming that the 
mass of the host star is equal to 0.93 Mq, the derived minimal 
mass of the inner planet is about of 0.065mj and its semi-major 
axis ^0.11 AU. Thus it would correspond to a hot Neptune- 
like planet. 
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Fig. 14. The elements of the best-fit solutions to the 2-planet model of the HD 188015 RV data projected onto the (P,e)- and 
(P, A')-plane. Upper plots are for the outer planet, bottom plots are for the inner companion. The (xl)^^^ of the best-fit solutions 
are marked by the size of symbols. Largest circle is for OtJ)'^^ = 1.021 (the best-fit solution, Table 1); smaller symbols are for 
Icr solutions with (xlY^^ e [1.02, 1.07), 2cr solutions with Oty)'''^ e [1.07, 1.14) and 3cr solutions with (^2)i/2 g [114^ ^ 25) 
(smallest, filled circles), respectively. 
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Fig. 19. Lomb-Scargle periodogram of the residuals of the 
long-period singal in the HD 1 14729 data. 
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Fig. 20. A histogram of (xlY^^ for Keplerian fits to 21,000 
sets of scrambled RV residuals of the synthetic signal of the 
outer planet in the HD 1 14729 planetary system (Table 1). The 
position of the best-fit to the real data is marked with an arrow. 



5. HD 147513 

The last exam ple which we exam ine is HD 147513. This star 
is observed bv lMavor et al.l ( l2004t) who published 30 measure- 
ments covering 1690 d. In this case, an rms of 5.7 m/s of the 
1 -planet fit is comparable with the mean of measurement er- 
rors, nevertheless our attention draws an irregular scatter of 
the data about the synthetic RV curve, Fig.l22h. The discovery 
team announced a Jupiter-like planet of minimal mass 1.21 mj 
in ^ 528 d orbit. The host star is very active (/?„„ = -4.4). 
Following Wright (20051) we estimate its jitter ~ 9 m/s. 



We searched for the second planet in short-period orbit, 
[2, 1 36] d. The results of the hybrid search obtained for 2-planet 
fits are illustrated in Fig. |^ for the statistics of the gathered 
best-fits and in Fig. I22h.b.c.d for the graphical illustration of 
the best-fit solution. These results remind us the HD 190360 
case. The (xl)^^^ of the best-fit has well defined minimum at 
f c - 7.894 d. Its rms ^4.1 m/s means a significant improve- 
ment of the single-planet fit. This is seen in Fig. 122b. d: appar- 
ently the synthetic RV curve is ideally close to the measure- 
ments. A relatively large ~ II m/s and mass of the host star 
1.11 Mq imply the minimal mass of the inner planet 0.1 15 iri} 
and semi-major axis about of 0.08 AU. 
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Fig. 17. The best-fit solutions to the HD 1 14729 RV data pubUshed in lButler et"aD(l2003l) . a) Synthetic signal of 1 -planet model, 
b) 2-Kepler model is consistent with the measurement error, (xlY^^ - 0.775, an rms ^ 3.6 m/s. c) Period-phased RV signal of 
the inner companion, d) A close-up of the synthetic signal of the best-fit 2-planet configuration (see Table 1). 
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Fig. 23. Lomb-Scargle periodogram of the residuals of the 
long-period signal visible in the RV data of HD 147513. 

The LS periodogram, Fig.|23j reveals features which recall 
us the first system examined in this paper, about HD 208487. 
The short period of ^ 7.9 d seems to be aliased with a much 
longer period of 810.8 d. Figure |23 is for the histogram of 
(xiy^^ computed for the scrambled residuals. The probability 
of a random character of the residual signal is negligible. At 
this time, such a characteristic is almost the same as for the 
HD 190360 system; in that case the 2-planet solution also pro- 
duces very small rms which is less than the joint measurement 
uncertainty. 

Still, mostly due to the activity of the star, the periodic- 
ity offlie_^Vvmations can be explained by jitter According 
to lMavoretal] ( l2004 . HD 147513 is a young G3/5V dwarf, 
only 0.3Gyr old. Its spectrum reveals a strong emission in the 
Ca II H line. Although the estimated short rotation period of 
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Fig. 24. A histogram of (xlY^^ for Keplerian fits to 30,000 
sets of scrambled RV residuals of the synthetic signal of the 
outer planet in the HD 147513 planetary system (Table 1). The 
position of the best-fit to the real data is marked with an arrow. 

the host star, about of 4.7 d jMavor et al.ll2004 . diff'ers from 
the secondary period visible in the RV signal (about of 7.9 d), 
it is suspiciously close to this approximation. However, such 
coinciden ce of periods still does not exclude a planetary ex- 
planation (iSantos et al. 2000, 2003) because the activity can be 
induced or amplified by relatively massive and close planetary 
companion. 

6. Conclusions 

It is extremely diflicult to interpret the RV measurements when 
their number is limited and their time coverage is poor We are 
looking for signals which have an amplitude comparable with 
noise. Although is akeady possible to pick up some promising 
extrasolar systems in which one planet is clearly visible and an 
another putative body is at the detection limit, a claim that the 
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Fig. 18. The parameters of the best-fits to the 2-planet model of the HD 114729 RV data projected onto the (P, e)- and {P, K)- 
plane. Upper plots are for the outer planet, bottom plots are for the inner companion. The values of (xl)^^^ of the best-fit solutions 
are marked by the size of symbols. Largest circle is for ixl)^^^ equal to 0.775 (the best-fit solution, Table 1); smaller symbols 
are for la solutions with {xiy"^ e [0.775,0.83), la solutions with (xlY'^ e [0.83,0.915) and 3(r solutions with (xiy^ e 
[0.915, 1.02) (smallest, fiUed circles), respectively. 
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Fig. 21. The projections of best-fit parameters to the 2-planet model of the HD 1475 1 3 RV data onto the {P, e)- and {P, /r)-plane. 
Upper plots are for the outer planet, bottom plots are for the inner companion. The values of (xl)^^^ of the best-fit solutions are 
marked by the size of symbols. Largest circle is for (xl)^^^ - 0.505 (the best-fit solution. Table 1); smaller symbols are for Icr 
solutions with (xiy^ e [0.505,0.619), 2cr solutions with 0^^2)1/2 g [0.62,0.773) and 3cr solutions with (xl)^'^ e [0.773,0.95) 
(smallest, filled circles), respectively. 
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for HD 147513. a) 1-planet model, 
parameters (K, P, e, o), - Tq) are: 



Fig. 22. The best-fit solutions to the RV data published in jMavor et all |2004 
The thick line is for the fit found with GA, assuming o-j - 9 m/s; its 

(28.40 m/s,520,49 d, 0.346, 288°.57, 3125.302 d), Tq^JD 2,450,000, 0^-2)1/2 ^ 9737^ rms=6.83 m/s. T he thi n 
line is for the fit without accounting for the jitter; its parameters are almost the same as in (Mavor et all I2OO4I) : 
(29.312 m/s, 528.1 12 d, 0.260, 282°.453, 2214.94 d), (xl)^^^ = 1.60, an rms = 7.26 m/s. b) 2-Kepler model is consistent with 
the measurement errors, Of^)'^^ - 0.505, an rms^ 4.1 m/s. c) Period-phased RV signal of the inner companion, d) A close-up of 
the synthetic signal of the best-fit 2-planet configuration. Table 1 . 



Table 1. Primary parameters of the model Eq. 1 and inferred elements of the best-fits found in this paper. To is JD 2,450,000. 
See the text and figures for error estimates. 





HD 188015 


HD 114729 


HD 147513 


HD 190360 


HD 208487 


Parameter 


c b 


c b 


c b 


c 


b 


c b 


M. [Mq] 


1.08 


0.93 


1.11 


0.95 




0.96 


m sin i [Mj] 


0.076 1.81 


0.066 0.88 


0.115 1.15 


0.064 


1.62 


0.121 0.46 


a [AU] 


0.08 1.22 


0.11 2.12 


0.08 1.35 


0.13 


3.93 


0.11 0.49 


Pld] 


8.389 476.23 


13.849 1171.19 


7.894 541.00 


17.104 2908.32 


14.496 129.51 


K [m/s] 


7.21 44.75 


5.87 17.71 


10.95 26.85 


5.20 


24.90 


10.47 19.08 


e 


0.23 0.24 


0.25 0.30 


0.30 0.09 


0.21 


0.37 


0.43 0.26 


CD [deg] 


103.9 242.4 


301.1 66.3 


9.1 316.8 


10.1 


10.8 


27.6 130.9 


Tp [JD-To] 


859.67 2270.10 


2496.77 1558.98 


1047.23 2248.44 


2370.38 3517.82 


2815.87 2957.74 




1.021 


0.775 


0.505 


0.680 




0.581 


D-j [m/s] 


4.0 


4.0 


9.0 


3.1 




3.0 


0=^ [m/s] 


4.1 


3.5 


4.4 


2.2 




5.4 


(o-j2 4-5=^)"- [m/s] 


5.7 


5.3 


10 


3.8 




6.2 


rms [m/s] 


4.93 


3.59 


4.07 


2.97 




2.82 


Vo [m/s] 


-17.45 


-3.22 


-4.85 


-4.81 




7.97 



residual signal has a planetary origin is always very risky. To 
confirm the existence of the putative hot-Ne ptune candidates, 
a deep astrophysical analysis is necessary dMavor & Oueloj 
ll99.5USantos et al.l2003U B utler et alJ20oi . by far beyond our 
technical abilities. That includes many additional spectroscopic 
and photometric observations. Nevertheless, the dynamical en- 
vironment of the putative new planetary companions, several 
recent discoveries of Neptune-like planets and cosmogonic the- 



ories supporting their origin, permit us to state and investigate 
a hypothesis of the planetary nature of the RV variability. 

Our results may be useful in some aspects. Even if the 
residual periodicities are of pure stellar origin, accounting them 
for may improve the single planet fits, because it makes it pos- 
sible to remove a contribution of not-random RV signal from 
the data. When we deal with a small number of measurements, 
the short-term contribution of jitter does not average out and 
it alters the long-period solution. Yet the stellar variations may 
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have origin in planetary-induced activi ty. This effect has been 
analyzed bv lSantos et alJ (|2p00, 2003). An elimination of the 
planetary hypothesis also requires additional observations and 
analysis. 

Basically, we use an independent, FFT-free method of the 
global analysis of the RV data and early detection of plausi- 
ble planets. It can be used for checking the results of 1-dim 
periodogram-based analysis, extending it into the multidimen- 
sional space of the orbital elements. If some of our results will 
be confirmed, a similar approach can be used to detect smaller 
bodies in short-period orbits in the other, already discovered 
extrasolar systems. Such studies can be useful for planning fu- 
ture observing sessions. 

We did a similar but very preliminary analysis of the RV ob- 
servations of a few other stars hosting single Jupiter-like plan- 
ets. For instance, we found a short-te rnivaria bility of the RV, 
about of 7 d in the HD 19994 system jMavor et al. 2004). The 
second companion would have K ~ 9 m/s, larger than the mean 
uncertainty of the measurements, about of 6.6 m/s. In the same 
paper, we found the RV data of HD 121504 and a curious RV 
variability which may be well modeled by a new 2-planet sys- 
tem close to the 6:1 mean motion resonance. Likely, the data 
collected and published by Doppler planet searches akeady 
hide many interesting and unusual planetary configurations. 
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